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Abstract 

A problem of amplitude reconstruction in terms of the given angu- 
lar distribution is considered. Solution of this problem is not unique. 
A class of amplitudes, correspondent to one and the same angular 
distribution, forms a region in projection onto a finite set of spherical 
harmonics. An explicit parametrization of a boundary of the region 
is obtained. A shape of the region of ambiguities is studied in partic- 
ular example. A scheme of partial-wave analysis, which describes all 
solutions in the limits of the region, is proposed. 

I. Let a quantum system, described by a wave function decay into 2 
spinless particles. We are studying angular distribution of decay products in 
the rest frame of initial state 

I(n) = |^(n)| 2 

and are going to obtain the wave function in the form of expansion by 
spherical harmonics 

*(™) = ^Ci m Yi m (n) . 

lm 

The solution of the problem: 



tt(n)| = Jl(H), *(n) = Jl{n)e^ H \ (1) 
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(f(n) is an arbitrary function of angles, 

c lm = j d 2 nY{ m (n)JT{Ky^ . (2) 

Substituting all possible functions <p(n) into this formula, one obtains the 
sets of coefficients Q m , correspondent to the same distribution I(n). In this 
the values of the coefficients span in the space {ci m } some region Q. Let us 
obtain a boundary of this region. 

We denote {lm) = a. Let us find the variation of coefficients when (p(n) 
changes 

On the boundary all 5c a are orthogonal to a single vector N a - normal to 
the boundary surface (fig.l). 

We introduce a scalar product in the space of coefficients in the following 

way: 

(a, b) = \ J2(a* a b a + a a b* a ) = E Re a a Re b a + E Im ajm b a . (3) 

a a a 

(N, 5c) = | E / d 2 n y/l(N^Y*e^ - N a Y a e^) 5^{n) = V 6<p(n) 
E KKe* - E N a Y a e-^ = e 2 ^ = S° 

a a l^a ly a 1 a\ 11 ') 



e i<p(n) _ _|_^ 



or 



-\E a N*Y*(n) 

(p(n) = axg^2 N a Y a (n) (+7r, for bottom sign). (4) 

This function y(n) is correspondent to the point, lying on the boundary of 
the region Q 

c p {N) = ±Jd 2 n ^I^Y^n)^^^^. (5) 

The sign ± in this formula is excessive. The region Q is symmetrical in 
reflections about the origin, because the angular distribution does not change 
in the replacement \& — > — \& (c — > — c). For the points on the boundary this 
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replacement is equivalent to the transformation iV — > —N. In fact, sign 
ambiguity in (|5]) corresponds to the replacement of outer normals to inner 
ones. 

Phase (£|) does not change in multiplication of all N a by positive number. 
One can choose the unit normal N a : (N, N) = \N a \ 2 — 1- 

a 

Let the unit normal N be given on a surface T. Thereby a mapping 
of the surface into unit sphere S is specified: 

r % s, n = N(c), c g r, n es . 

This mapping is called Gaussian mapping [Q] . It possesses the property 
equivalent to the definition: 

(N, dc) = 0. 

Formula (H) specifies a parametrization of the surface T = dQ, the 
parameter is unit normal N. Thus T is defined by the inverse Gaussian 
mapping 

S^T, c = c(N), N G S, c G r . 

In order to obtain a projection of the region Q in some subspace {c/ m , I = 
O..L max }, one should restrict N a in formula (Q) into this subspace. Values c@ 
obtained specify a boundary of the projection of the region (fig. 2). 

In order to obtain a slice of the region Q by specified hyperplane {ci m = 
c^l - fixed numbers , / > L max }, one should solve a system of equations 
on N Q 

[d 2 nVlY*e ia *^ N ° Y <* = cf. 
This problem is more difficult than preceding one. 



Example. Let 1=4-. What initial states (except the obvious S-wave 
\& = ^7= ) are able to give this distribution? 

The answer: \I/ = -^^^ '■ The boundary of the region f2 is given by the 
formula 

C/3 = — L J d 2 n Yje ia **'E NaYa . (6) 
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Now the problem consists in clear representation of this region. 

1) Projection of Q in finite dimensional space {ci m , I = O..L max } lies inside 
unit sphere 

^max CO p 

E C « C «<E C « C « = d 2 n\V\ 2 = l. 
o o J 

2) Projection of Q on each wave LM covers a circle inside unit circle on 
Argand plot: 



when Ni 
1 



e^, I = L, m = M 
■ a \ ~ . ar 



clm 



0, otherwise 
/ d 2 n |y iM |e 



» arg Y LM _ & i arg F^m 



g E N a Y a = 4> + arg Ylm, 

1 



47T 



/ d 2 n \Y LM \ 



(L,M) : 


(0,0) 


(1,0) 


(1,1) 


(2,0) 


(2,1) 


(2,2) 


= \e** J dcos9 ■ 
-l 

= e i(t> - 


1 

i; 
1-; 


\/3| cos 

2 ' 

0.962, 


y/f sin6> 

4 V 2' 

0.866; 


^|3cos 2 6>-l| 

2 /6 

3 V 3' 

0.861, 


y/^| sin6> cos 6* 
/? 

V 6' 

0.913, 


S m 2 6 
fi 

V 6 

0.913 



(cl -m = cla/) 
3) Projection on S,P -waves (a space {c ,c w }). 

When rotating phases of all N a by O , the phases c a also rotate by 4>o- 



N n -> N n e i<l>0 



i<k> 



In multiplication of iV^ by a positive number the c a do not change. 

Let Nq = 1. We will measure the phase of P-wave with respect to the 
phase of S-wave 

c 10 ^c 10 e- iar s co ,c ^ |c |. (7) 
Let us denote v^^io = P e ^ 

i arg(l+pe l< ^ cos 9) 



1 + p COS I 



a/1 + 2p cos 9 cos 4> + p 2 cos 2 



From 



1 /" 27r , f 1 + sce*^ 

c = — / dip ■ 

47T J0 J-p p \/l + 2xcos0 + x 2 



x = p cos 6* 
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sin /-(p+cos 4>)/ sin ^ i _ e i<P cos + e ^ sin • m x + cos . 

aw ; u = — - — — 



2p J(-p+cos0)/sin^ | sin0| Vl + u 2 



sin < 



e %4> ^ 
2^ 



\^Jl + 2pcos0 + p 2 - y^l - 2pcos0 + p 2 (8) 



.. . p + cos + e VI + 2p cos + p 2 \ . . 

-2 sin In ■ — ■======) e=sgnsm0. 

— p + cos + e VI — 2p cos + p 2 / 



Analogously 



V3 f 2w fP dx 1 + xe^ 

c io = -r- / "V 5 ' / -5- z 



47T Jo i-p p 2 VI + 2XCOS0 + X 2 



V3 



■-,2 



2e^((p + cos 0)^1 + 2pcos0 + p 2 



— (— p + COS (j))\Jl — 2p COS + p 2 ) 

-4e 2i *(^/l + 2p cos + p 2 - ^1 -2pcos0 + p 2 ) (9) 

p + cos + eVl + 2pcos0 + p 2 \ 



+i(3e^ + l)|sin0|ln 



— p + cos + eVl — 2pcos0 + p 2 



The limits 



p — > : Co — > 1, Cio — > (normal lies in the plane of S-wave) 

p — > 00 : c — > 0, Cio — > ^e 1 ^ (normal lies in the plane of P-wave) 

Fig. 3a shows the projection of the region Q into coordinates (c , | c 10 j ) - 
Dotted band is the projection of the boundary <9f2, the projection of Q is 
shown by sloppy hatching. Fig. 3b shows the slices of dfl by level surfaces 
Co = Const in projection on Argand plot (Re cio, Im cio). When Co — > 1, the 
slices contract into a point (0, 0) (the P-wave contribution disappears). When 
c — > 0, the slice reaches the circle with the radius a/3/2. For intermediate 
values < c < 1 the slices are of the oval shape. In this the length of 
radius- vector |cio| changes from minimal value to maximal value in the band, 
displayed on fig. 3a. The slices of Q by surfaces Co = Const lie inside the ovals. 
Fig. 3c shows 3-dimensional region of ambiguity. 
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Expansions of the functions \l/(n) = e %v ^ n > /a/47t in common contain in- 
finite number of harmonics. Figures display the projections (shades), which Q 
casts from infinite dimensional coefficient space in specified 2- or 3-dimensional 
subspaces. A contribution of the infinite number of harmonics not involved 
in these figures can be evaluated. We use the equality 

oo „ 

EM 2 = d2fi \^\ 2 = 1 E |c a | 2 = i-c 2 -|ci | 2 . 

On fig. 3a two circles with radii 0.95 and 1 are shown. For the points of 
Q, lying between these two circles, the contribution of "invisible" harmonics 
does not exceed 10% of total probability 

0.9 < c 2 + |ci | 2 < 1 M 2 <0.1 

a^S,P 

The same region is hatched on fig.3b. On the bottom figure this region is 
correspondent to a part of Q, lying outside a sphere with radius 0.95. 

The figures are obtained as follows. 10 6 random points (p, <ft) were 
generated with uniform distribution in rectangle 

{0 < p < 30, < <p < 2vr}. 

Coefficients Co , cio were calculated by the expressions @ , (|9|) . Redefi- 
nition (Q) was performed: cio — > cioCq/|co|, cq — > |co|. Values (cq, \c\q\) 
were placed in 2-dimensional histogram. Besides, 5 distributions were 
filled on Argand plot (Re cio, Im cio) for Co, lying in the intervals 

{p - 0.03 < c < p + 0.03} p = 0.1, 0.3, 0.5, 0.7, 0.9 

4) One more symmetry of the region Q. 

In transformations of Ni m by matrices T> l mm , (Wigner's functions) the 
coefficients Q m are transformed by the same matrices 

Ni m ->• ^V l mm ,(u))Ni m , => ci m ^^2V l mm ,(u)ci m i 

m' m! 

(proved by the replacement of integration parameters in (|6|)). 

Unitary matrices T> l mm , conserve scalar product (|3|). In T> l mm , -rot at ions 
the region Q transforms to itself, with its normals. 
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If one has studied in detail the region Q in the considered example, then 
the regions of solutions for non-isotropic distributions / 7^ Const can also be 
obtained. Formula (||) is linear to \fl: 

if Vl = cti\fh + at2\fl2, then cp = a\cp + c^c^. 

Expanding \fl by spherical harmonics: \fl = b a Y a and using in (^) the 
multiplication theorem for spherical functions 

Y a Yp = linear combination of Y2 , 

we see, that the region fli^const can be obtained from the region fli=c nst by 
a linear transformation with coefficients depending on the distribution I. 

Substantial questions 

1. For the scalar product (Q) we have chosen the Euclidean scalar product 
in real space, formed by real and imaginary parts of c a . (This scalar product 
is also equal to real part of the Hermitean scalar product in complex space 
c a .) It is clear that the scalar product is an additional structure, which 
does not influence the shape of the region Q. One can easily show that 
the replacement of scalar product leads to some linear transformation of the 
parameter N — > N(N) in (||), i.e. the shape of the region conserves, only its 
parametrization changes. One can imagine this in the following way. The 
replacement of the scalar product changes the orthogonality definition and, 
consequently, the definition of the normal to the surface. When the scalar 
product replaces, the field of normals N on the boundary dQ changes, but 
the surface dfl remains the same. 

2. In infinite dimensional space of coefficients c a the set Q is not a region. 

00 

The f2 is a surface, lying on a sphere J2 \ c a\ 2 — J d 2 n I(n) = 1. This surface 



maps to the region in projections into any less space. Expressions (0) and 

(|5[) in initial space coincide (Q = dQ), because any function <p(n) (for which 

• /-\ ' 00 

e vp{n) j s quadratically integrable) can be presented in the form arg N a Y a (n). 



Actually, we consider the projection of Q in some finite dimensional space 
from very beginning. In derivation of formula (|5p the vector N a was restricted 
to this space. 
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3. In those points of a sphere no, where Y^N a Y a (fio) = 0, the function 
<p*(n) = argj] N a Y a (n) has a break (branching point). The correspondent 
wave function ^ = \fle lLp * also breaks in this pointQ, except a special case 
/(no) = 0. This circumstance is complicated by the fact, that all functions 
of form (P, close to Nl/*, also have breaks. (The phase of ^* changes by 27r in 
bypassing around n . All functions, close to ^P*, also possess this property. 
If l^*] 2 = l^] 2 = / 7^ in the vicinity of no, then \I/ should break in this 
vicinity.) 

Therefore, Q includes regions, correspondent to irregular functions. If 
the regular functions $ should be found, for which |\I/| 2 = / exactly, then 
additional analysis is needed to discard the regions of irregularity. Outside 
class ([I]) smooth functions \l/ exist, which are close to ^ everywhere, except 
a small vicinity of n : 

Ve,5 3f GC 1 ^ 2 ) |*(n)-^*(n)| < e Vn : |n-n | > 6. 

Inside 5-vicinity of n the difference ^— ^* is bounded. (We present functions 
^ = z/\z\ and \I/ = z/(|z| + 5) as an example, see fig. 4). 

The phase of \I/ changes by 27r in bypassing around n - From the conti- 
nuity, a point exists in 5- vicinity of n , where the \l/ and the correspondent 
density |\1/| 2 vanish. Outside 5-vicinity |\l/| 2 can be close to I as one likes. 
Coefficients c a for \l/ are close to the coefficients for ty*. Therefore, the dis- 
tinction of and is not physically significant. 

II. Let us consider a fermion, decaying into 2 particles with spins and 
|. The initial state is described by 2-component spinor 

The angular distribution observed is 

7(n) = |^i(n)| 2 + \^ 2 (n)\ 2 - 

The space of states is a direct product of an orbital space of decay 
particles and a spin space of the final particle with the spin |. 

however, it is single- valued, bounded and expansible by spherical harmonics in the 
vicinity of no- A detailed study shows that formula rtq) holds true for these functions. 



8 



The state with definite orbital moment I and projections m, X is: 

\lm; \X^ = Y lm (n)x\, X + i = f q ) ! X -\ = ( 1 ) ' 
The state with definite total moment J and projection M is: 

\ JM ) = J2 (l M — X; \X\JM^j ■ \l M — X; ^Xj. 
i=J±\ 

If the coefficients of expansion of each component of the spinor \& by 
spherical harmonics are known 

ImX 

it is possible to obtain spin composition of the initial state 



CJM 



Let us find ^f(n) 



JM\*)= Yl cl M _ x (lM-X;\X\JM 



i=J±\ 

A=±i 



I ^1(^)1 = \J cos0(n), |^2(^)| = \Jl{n) sin0(n), 

0(n) is an arbitrary function of angles. This function defines a degree of 
polarization of the final particle with the spin |, which is not measured in 
the experiment. 



mJn) = J 1(H) cos <p(n)e iipl{H \ * 2 (n) = Jl(n) sin 0(n)e^ 2 ^ 



d*nY L (n)^) e ^-{2f$; \Z\ (10) 

One can apply to this formula the analysis described above. The result is 
the following. The boundary of the region of ambiguity is parametrized by 
complex numbers N l a . The point on the boundary c l a (N) is defined by the 
formula (|10|) , where 



<fi = arg Nl,Y a and tan 



V N 2 Y I 
v N l Y I 
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III. Let the process be described by density matrix. We suppose that only 
spinor components were mixed, i.e. only the information about polarizations 
is lost. For example, consider the process: + | — > X — > + |, 2 plane 
waves scatter in the initial state, the initial fermion is not polarized, the 
polarization of final fermion is not measured. 

The initial state is 



= 5{k-kiz)x\i ~ E \ -7~ Yl0 ^ ' *V 



i 



47T 

The orbital part of the wave function \ij is known, the spin part is unknown. 
The initial state is described by density matrix: 



Ai=±2 



The intermediate state X depends on A^: I^a,/ = M\i\ij, where M is a 
reaction amplitude. Components of the wave function X and coefficients of 
their expansions by spherical harmonics represent matrix elements of reaction 
amplitude 



^XiX f (n) = (n\ f \M\i\i), \n X f j = 5 (ft - n )xx f , 
- ^lm; \\f\M\i\i\ \lm;\\ f \ = Y lm (n)x\ r 



Hm 

The angular distribution observed is 



I (ft) = J2 (^A/| MpiM + \n\j 
1 s " ' 

A/=±2 Px 



= \ E {^f\M\i\) {iMM + \n\ f ) = \ £ I^A^fOI 2 

A/=±| X i =± \ 

AlA/V ; V v; (sm^i(n), \ = -\ \sm0 2 (n), X f = - 

c X a iXf = Jd 2 nY:(n)^ XiXf ( 



f (n). 
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This formula contains 6 arbitrary functions. For the points on the boundary 
these functions have a form: 



<p XiXf = arg£ N^Y a 

a 

R-+ 2 + R-~ 2 -R + - 2 ~R ++2 ±^D 
2(R-+R++ + R—R+-) 

-r-+ 2 + r— 2 + r + - 2 -r ++2 ±Vd 



tan 



tan02 2(Rr+R- + R+-R++) 

D = (R ++2 + R+- 2 -R— 2 -R- +2 ) 2 + 4(R- + R ++ + R—R+-) 2 , 

where R XiX f = \J2 N a XfY a\- 

a 

Notions. 



1) In view of axial symmetry of the problem the distribution I does 
not depend on azimuthal angle. The amplitude \& is also proposed to be 
axially symmetrical, its expansion contains harmonics Yi only. Derivation 
of ambiguity region in the space {q } is equivalent to sectioning of fl by this 
space. The result is described by the same expressions, in which the indices 
a are restricted to the set (70). 

2) All mentioned above corresponds to the case of inelastic scattering 
(particles in initial and final states are different). For elastic scattering it is 
necessary to extract explicitly in S-matrix a contribution of the non-scattered 
wave 

S = 1 + iM. 

The consequence of unitarity of S'-matrix is an optical theorem for the am- 
plitude of elastic scattering: 

Im*U =i ~ o-tot- 

(If the initial state is fixed, then the optical theorem is the unique consequence 
of the unitarity of S.) This condition define a hyperplane in the space of 
coefficients. The slice of the region of ambiguity by this hyperplane is the 
solution of the problem in the case of elastic scattering. 

3) For pure elastic scattering the combined imposition of unitarity and 
rotational invariance of S'-matrix fixes the continuous ambiguity (for example, 
see 0] §125). With nonelastic channels this does not occur. Rotational 
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invariance leads to diagonality of S'-matrix with respect to quantum numbers 
of the total moment 

(j'M'n'\S\JMn) = 8jj.8 MM .S^ , 

index n includes all other quantum numbers, e.g. the type of a particle. 
The matrices Sj^Jf^ are unitary. This condition provides non-linear rela- 
tions among partial amplitudes of different reactions. An exact inclusion of 
the unitarity condition is possible only in the measurement and combined 
analysis of the whole variety of connected channels. 



An extensive literature is devoted to the question of continuous ambigui- 
ties in phase analysis (see review in ||, pp. 291-308). The reaction amplitude 
is ambiguously reconstructed from experimentally observable quantities, be- 
cause the part of information about amplitude is lost in the measurement. 
Non-observable quantities (phases, polarizations) can be functions of angles, 
energy, etc., therefore the final expression for the amplitude contains func- 
tional arbitrariness. This arbitrariness is significant (if one substitutes fast 
changing phase (p(n) in formula fl5|), then the contribution of low harmonics 
decreases, the same distribution is described by high harmonics). One can 
diminish the ambiguity only imposing extra restrictions on the form of ampli- 
tude, i.e. adopting the particular model. Note, that in the case of the known 
region of ambiguity the inclusion of extra restrictions reduces to sectioning 
of the region by the surface of the restriction. 

For example, it is known that the initial state has definite parity. 
In this case one should consider the region Q in the space / have 

the given parity }. In this one can use formula (|5|), restricting I in the 
indices a and (3 to the set with the given parity. One can prove this 
considering the slice of ft, by this space or repeating the reasonings in 
section I for Vl/ with the definite parity and restricting the integrations 
on semi-spheres. 

The state ty(n) = ^Jl(n)e vp ^ includes in general an infinite number of har- 
monics. Coefficients q m tend to zero at I — > oo (fastly, if ty(n) is sufficiently 
smooth). In practical partial wave analysis the experimental distribution is 

fitted by the finite segment of the expansion: ^ = J2 ci m Yi m , in this the 

o 
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finite number of solutions has been obtained J3|. This happens because the 
segment of the expansion approximates the exact solution, and just for the 
finite number of values of parameters this approximation is the best (fig. 5). 

At L max — » oo we turn back in infinite dimensional space, where the 
number of solutions is infinite. Apparently, when L max increases, the fit 
becomes unstable or new solutions appear, which fill the region Q in the 
limit L max -> oo. 

There is much speculation that the consideration of high harmonics is 
not necessary, because noise exceeds a signal in high harmonics due to finite 
angular resolution of the equipment and limited statistics. The contributions 
of high harmonics cannot be precisely measured, but certainly this does not 
mean that they are equal to zero. We have observed that weak and sharp 
conditions 

cj m — ► / — > oo and Q m = I > L max 

have principal difference. In the first case the large "ambiguity island" ex- 
istsQ. In the second case the finite number of solutions will be found, closest 
to the boundary of the island. The situations are possible when resonant 
structures in low harmonics can be also described by correlated contribu- 
tions of high harmonics. One can avoid this uncertainity, if the whole region 
of ambiguity will be found in analysis. The general scheme of such analysis 
is the following: 

• Energy axis of the system considered is divided in bins, the angular dis- 
tributions are filled in each bin. These distributions are approximated 
by smooth functions, e.g. spherical harmonics: 

J (X> = XXm^m(w), ti m = J d 2 n I(n)Y lm (n) = (Y lm 

lm 

• The boundary of the ambiguity region is obtained by formula (H). The 
integral on a sphere entering this formula can be numerically evaluated. 

This integral is not expectation value of some quantity, in its calculation the 
whole distribution lift] must be known. Also, it should not be calculated 



taking a root y/I(n) from each bin content and summing over bins. Many- 
dimensional distributions are rarely filled and the errors of such integration 



2 The region of ambiguity presented on fig. 3 fills almost the whole unit sphere. 
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will be great. The approximation of I(n) by the smooth function in the 
previous step is necessary. 

A many-dimensional surface obtained can be kept in the form of linear 
frame. A displaying technique for many- dimensional surfaces presented 
in this form is being developed nowdays ||. 

Note, that our approach does not require the fitting of the experimental dis- 
tribution. This excludes various problems, such as the derivation of all best 
approximations, which requires either repeated applications of fitting proce- 
dure with different start points or complex algebra for discovering analytic 
relations between different solutions. 

A presence of resonance in some wave is established independently on a 
model only if the projection of ambiguity region in 3-dimensional space 

Argand plot for this wave x energy axis 

is a narrow tube containing only resonant paths (when the energy changes, 
the absolute value of the wave has a bump, the phase changes by n). In 
other cases the non-resonant paths are present along with resonant ones. In 
this case one may only state that the experimental data do not contradict 
the presence of the resonance in the given wave. 

A model predicts the path c(E) (or the class of paths). If the ambiguity 
region is obtained, one can easily check, whether the path is contained in the 
ambiguity region and also find the area of values of model parameters, when 
the paths are contained in Q(E). Particularly, these parameters might be 
masses and widths of resonant states. 



14 



References 

[1] Francis G.K. A Topological Picturebook. Springer- Verlag 1987,1988; 
(Moscow, Mir, 1991, Russian). 

[2] Landau L.D., Lifshits E.M. Theoretical Physics V.3: Quantum Mechan- 
ics. Moscow, Nauka, 1974. 

[3] Nikitiu F. Phase Analysis in the Physics of Nuclear Interactions. 
Moscow, Mir, 1983 (Russian). 

[4] Sadovsky S.A. Preprint IHEP 91-75, 1991; Prokoshkin Yu.D., Sadovsky 
S.A. Preprint IHEP 92-16, 1992 (Russian). 

[5] Nikitin I.N., Talanov V.V. Preprint IHEP 94-13, 1994. 

Received December 8, 1994- 



15 



Figures captions 

Fig.l. For the point inside the region Q, the linear span of all possible vari- 
ations 5c a coincides with the whole space. For the point on the boundary of 
Q all Sc a lie in the tangent plane to the boundary surface and are orthogonal 
to normal vector N a . 

Fig. 2. The region f2, shade and slice. 

Fig. 3. The ambiguity region for isotropic distribution in projection on 
SPo- waves (see text for explanations). 

Fig. 4. Function ty* has a break in the point Hq. All functions, close to fy*, 
for which |\&| 2 = |\I/*| 2 , have breaks in the vicinity of Hq. Continuous functions 
\l/ exist, which are close to ^ everywhere, except the small vicinity of n - 
For these functions the |\l/| 2 vanishes in some point inside the vicinity. 

Fig. 5. The restriction of the number of harmonics and model substitutions 
for the amplitude lead to the fact that the solution is considered in some 
narrow class of functions. This class and the set of exact solutions Q in 
general do not intersect each other. The distance between points in these 
two classes become minimal in a discrete set of points. 



16 




fig.l 





fig. 





fig. 4 




fig. 5 



